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LONG-TERM  GOALS 

The  long-term  goal  of  this  effort  is  to  demonstrate  the  applicability  and  performance  of  the  Delft3D 
hydrodynamic/morphodynamic  software  system  to  the  case  of  the  hydrodynamics  and  contaminants 
dispersal  and  advection.  Delft3D  has  until  now  been  applied  extensively  to  strongly  astronomical, 
wind  and  wave-forced  (shelf  sea  and  nearshore)  conditions  but  not  to  the  weakly-forced  conditions 
such  as  occur  off  the  California  coast.  This  is  of  interest  to  the  Navy  because  of  naval  operations  and 
the  presence  of  a  Navy  SEAL  training  ground  just  north  of  the  Tijuana  River  outlet. 

In  a  separate  task,  the  performance  of  Deltares’  river  model  SOBER  is  evaluated  for  the  case  of  the 
Pearl  River  (MS, LA),  also  of  interest  to  the  Navy  with  respect  to  training  grounds  near  Stennis  Space 
Center  (MS). 

OBJECTIVE 

The  objective  of  this  study  is  to  demonstrate  the  capabilities  of  the  Delft3D  hydrodynamic  system  by 
developing  a  multi-domain  (nested)  cascade  of  models  ranging  from  the  Southern  California  Bight 
scale  to  the  nearshore  area  around  the  Tijuana  River  mouth  and  Imperial  Beach.  The  model  is  to  be 
driven  with  boundary  conditions  from  a  larger  ocean  shelf  model  (wind,  tides,  salinity,  temperature, 
etc.)  for  a  number  of  rainfall  events  (scale  of  days  and  weeks),  and  to  be  validated  against  measured 
data  inside  the  domain. 

An  analysis  will  be  made  of  what  conditions  determine  the  trapping  and  mixing  of  the  plume  inside  the 
surfzone  and/or  the  propagation  of  the  plume  through  the  breakers  and  onto  the  coastal  shelf,  and  - 
through  a  data  denial  study  -  what  physical  processes  are  important. 

For  the  Pearl  River,  a  demonstration  will  be  given  of  the  performance  of  the  SOBER  river  model. 
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APPROACH 


The  work  is  divided  into  4  tasks 

1.  Application  of  the  Delft3D  model  to  the  San  Diego  Bay  and  coastal  region 

2.  Enhancement  of  the  Delft3D  system  for  operational  use 

3.  Peal  River  modeling  using  SOBER 

4.  Development  of  beta-version  of  San  Diego  Delft3D  model  in  nowcast/forecast  model 

WORK  COMPLETED 

1 .  Application  of  the  Delft3D  model  to  the  San  Diego  Bay  and  coastal  region 

We  have  set  up  a  Delft3D  sea  shelf  model,  coupled  it  with  a  larger-scale  model  and/or  data  to  provide 
boundary  conditions,  including  flow  boundary  conditions  and  solar  radiation.  In  this  way  the 
hydrodynamics  of  the  eddies  in  the  San  Diego  Bight  can  be  simulated.  We  use  point-source 
contaminants  at  the  Tijuana  River  Mouth  as  input  in  order  to  perform  simulations  of  the  plume 
contaminant  dispersal  and  water  quality  conditions  on  the  San  Diego  Bight.  Specifically,  we  performed 
1 -month  model  simulations  at  the  time  of  rainfall  events  (February  2009)  to  analyze  model  success  at 
simulating  local  processes. 

Only  recently  COAMPS  and  NCOM  data  has  become  available  for  this  period  which  will  be  used  to 
drive  the  already  set  up  model.  We  will  show  the  deterioration  of  the  performance  of  the  model 
through  a  data-denial  study,  will  be  analyzed  in  order  to  determine  the  optimal  input  of  data.  Impacts 
of  mixing  of  the  contaminants  across  the  surf  zone  will  be  examined.  It  is  expected  that  in  high-surf, 
in  combination  with  the  low  river  jet  momentum  conditions,  the  plume  will  be  trapped  within  the  surf. 
In  contrast,  it  is  expected  that  in  higher  river  jet  momentum  conditions,  the  plume  will  quickly  cross 
the  surf.  Additional  complexity  can  occur  when  high  shear  is  present  between  the  wave-driven  long¬ 
shore  currents  within  the  surf  and  offshore  wind/tidal/pressure  gradient  forced  currents  move  in 
opposite  directions. Scientific  or  technological  objectives  of  this  effort. 

In  order  to  do  the  analysis  a  “drifter”  code  modification  has  been  made  in  the  code  which  allows 
numerical  drifters  to  move  horizontally  in  water  column  along  a  prescribed  path  in  order  to  simulate 
REMUS  trajectories.  At  each  horizontal  position  all  variables  in  the  entire  vertical  will  be  stored. 

Given  a  known  vertical  position  of  the  instrument  the  local  value  of  any  variable  can  be  retrieved  from 
this  column  of  data 

2.  Enhancement  of  the  Delft3D  system  for  operational  use 

We  have  implemented  a  roller-model  routine  in  Delft3D-FLOW,  which  is  able  to  calculate  the  wave 
directions  independently  of  SWAN  on  the  FLOW  grid.  This  module  has  been  tested  extensively. 

This  calendar  year  a  mass  balance  diagnostics  module  will  be  added  to  the  Delft3D  code. 

3.  Pearl  River  modeling  using  SOBER 

Three  visits  have  been  made  to  NRL-Stennis  to  demonstrate  the  capabilities  of  the  SOBER  model,  and 
to  make  an  implementation  of  the  river  shed  using  available  online  data  and  discharge  events. 
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The  model  set-up  is  conducted  with  the  two  hydrodynamic  modules  “Flow”  (ID)  and  “Overland 
Flow”  (2D).  For  demonstrating  the  range  of  possible  modelling  options  in  the  software,  we  set-up 
different  models:  i)  ID  only  model,  ii)  2D  only  model,  iii)  a  combined  1D2D  model.  Results  are  still 
forthcoming. 

4.  Development  of  beta-version  of  San  Diego  Delft3D  model  in  nowcast/forecast  mode 
This  work  is  projected  for  the  coming  year. 

RESULTS 

The  curvilinear  mesh  of  the  San  Diego  Bight  covers  the  30  km  stretch  of  coast  from  Rosarito  (Mexico) 
to  Solana  Beach  see  Figure  1.  In  order  to  capture  the  temperature  and  salinity  stratification  in  San 
Diego  Bight,  it  is  necessary  to  run  the  model  in  3D  mode.  For  this  purpose,  20  vertical  sigma  layers  are 
defined.  The  largest  vertical  density  gradients  in  the  model  are  expected  to  occur  close  to  the  water 
surface  and  this  is  where  the  highest  vertical  resolution  is  required. 

Setting  up  boundary  conditions  for  the  model  has  turned  out  to  be  one  of  the  most  challenging 
problems  in  this  project  so  far.  Multiple  configurations  of  boundary  condition  types,  such  as  water 
levels,  velocity  and  Neumann  (water  level  gradient)  forcing  have  been  tested  along  the  three  (southern, 
western  and  northern)  open  boundaries.  The  application  of  water  level  boundaries  in  the  San  Diego 
Bight  model  often  gives  rise  to  large  instabilities  in  the  form  of  high  velocity  water  jets  entering  the 
model  domain.  It  is  thought  that  the  strong  vertical  density  variation  in  combination  with  relatively 
weak  forcing  agents  (weak  tidal  currents,  low  wind  velocities)  causes  these  instabilities.  The  jets  cause 
unrealistic  vertical  mixing  and  large  velocities  inside  the  model  domain.  Using  current  and  Neumann 
boundaries  conditions  along  the  lateral  boundaries  generally  causes  unrealistic  upwelling  close  to  the 
shore. 

An  important  finding  is  that  only  the  use  of  Riemann  (weakly-reflective)  boundaries  yielded  both 
stable  and  realistic  results.  Although  these  are  more  cumbersome  to  set  up  (one  needs  both  normal 
velocities  and  water  levels  at  each  boundary)  and  to  some  extent  cause  some  loss  of  control  over  the 
water  level  inside  the  model  domain. 
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Figure  1.  Computational  grid  San  Diego  Bight 

It  is  important  to  get  an  estimate  of  what  the  model’s  sensitivity  to  different  input  parameters.  For  this 
purpose,  a  sensitivity  analysis  has  been  carried  out  in  which  the  following  changes  have  been  made 
with  respect  to  a  reference  simulation 

•  No  wind 

•  No  heat  exchange  with  the  atmosphere 

•  No  currents  on  the  southern  boundary 

Results  are  shown  in  the  following  figure  and  indicate  that  it  is  important  to  have  wind  forcing  ,  heat 
exchange  and  currents  on  the  boundary  all  turned  on  in  the  model,  i.e.  no  process  is  dominant. 
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IMPACT/APPLICATIONS 

The  model  results  so  far  show  that  the  Delft3D  model  is  capable  of  simulation  the  essential  processes 
in  the  San  Diego  Bight  and  can  be  forced  by  an  outer  model,  and  wind  forcing.  This  tool  has  the 
potential  to  deliver  real  time  forecasting/nowcasting  capabilities  in  this  region. 

TRANSITIONS 

The  Delft3D  model  software  and  code  is  accessible  to  ONR  under  the  existing  DELFT3D  license  of 
contract  N00014-02-C-075,  dated  October  3,  2002. 

RELATED  PROJECTS 
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